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(57) Abstract 

A magnetic resonance (MR) imaging ap- 
paratus and technique exploits spatial infor- 
mation inherent in a surface coil array to in- 
crease MR image acquisition speed, resolution 
and/or field of view. Partial signals are bc- 
quircd simultaneously in the component coils 
of the array and formed into two or more sig- 
nals conesponding to orthogonal spatial rep- 
resentations. In a Fourier embodunent, lines 
of the *-space matrix required for image pro- 
duction are formed using a set of separate, 
preferably linear combinations of the compo- 
nent coil signals to substitute for spatial mod- 
ulations noimally produced by phase encoding 
gradients. The signal combining may proceed 
in a parallel or flow-through fashion, or as 
post-processing* which in either case reduces 
the need for time-consuming gradient switch- 
ing and expensive fast magnet arrangements. In 
the post-processing approach, stored signals are 
combined after the fact to yield the full data ma- 
trix. In the flow-through approach, a plug-in 
unit consisting of a coil array with an on boanl 
processor outputs two or more sets of combined 
spatial signals for each spin conditioning cy- 
cle, each directly corresponding to a distinct 
line in k-space. This partially parallel imag- 

ii^toll toe wid, a protows thr««oil array, and lanscr fiictore a« achievable with other co.l anangemenU. 
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Simultane usAcquisiti n of Spatial Harmonics (SMASH): 
Ultra-Fast Imaging with Radiofrequency Coil Arrays 



Background 

5 This invention is generally related to magnetic resonance imaging (MRI) using 

nuclear magnetic resonance (NMR) phenomena It is particularly directed to a method 
and corresponding apparatus for more efficiently capturing and providing MRI data 
suitable for use in multi-dimensional imaging processes. 

MRI is by now a widely accepted, medically imiwrtant and commercially viable 

1 0 technique for obtaining digitized video images representative of internal body tissue and 
structures. There are many commercially available approaches and there have been 
numerous publications describing these and other approaches to MRI. Many of these 
use multi-dimensional Fourier transformation techniques which are now well-known to 
those skilled in this art 

15 In general, MRI devices establish a constant homogeneous magnetic field, apply 

a specific additional bias field gradient in a known plane or region under consideration 
to orient nuclear spins, and apply a radiofrequency pulse or a sequence of pulses to 
perturb the nuclei. These nuclei in the known bias field gradient emit an RF signal in a 
specific band determined by the magnetic field distribution, and these RF emissions are 

20 detected by receiving coils and stored as a line of information in a data matrix known as 
the k-space matrix. The full matrix is buih up by successive cycles of conditioning the 
spins, perturbing them, and collecting RF emissions. An image is then generated fiom 
this matrix by Fourier transformation, which converts the fiequency information present 
in the RF oscillations to spatial information representing the distribution of nuclear spins 

25 in tissue or other imaged material. 

Magnetic resonance imaging has proven to be a valuable clinical diagnostic tool 
in a wide range of organ systems and pathophysiologic processes. Botii anatomic and 
fimctional information can be gleaned from tiie MR data, and new appUcations continue 
to develop witii each improvement in basic imaging technique and technology. As 

30 technologic advances have improved achievable spatial resolution, for example, 
increasingly finer anatomic details have been able to be imaged and evaluated using MR. 
At tfie same time, fast imaging sequences have reduced imaging times to such an extent 
tiiat many moving structures can now be visualized without significant motion artifacts. 
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Often, however, there is a tradeoff between spati^ resolution and imaging time, 
since higher resolution images require a longer acquisition time. This balance between 
spatial and temporal resolution is particularly unportant in cardiac MR, where fine 
details of coronary artery anatomy, for example, must be discerned on the surface of a 
5 rapidly beating heart. A high-resolution image acquired over a large fraction of the 
cardiac cycle will be blurred and distorted by bulk cardiac motion, whereas a very fast 
image may not have the resolution necessary to trace the course and patency of coronary 
arteries. Some of the fastest imaging sequences cunrently implemented, such as echo 
planar imaging (EPI), approach the goal of yielding images of reasonable resolution in a 

10 suitably short fraction of the cardiac cycle. Other approaches have also been tried to 
eluninate the effects of cardiac motion, including A-space segmentation, m which image 
acquisition is divided up over several cardiac cycles with ECG gating to ensure that the 
heart is in the same phase of systole or diastole during acquisition of each segment. 
Cine fan^es of multiple cardiac phases may be pieced together with this technique, with 

15 partial acquisitions for different phases occurring in each cardiac cycle. One problem 
with this class of techniques is that respiratory motion can change the position of the 
heart over the course of several cardiac cycles. Partial acquisitions will then be 
misregistered, and artifacts will result. In an attempt to eliminate or adjust for 
respiratory motion, breath holds, respiratory gating, and navigator echo gating 

20 techniques have all been tried, and each of these techniques has had some significant 
successes. Nevertheless, an imaging strategy v*ich allowed high-resolution images to 
be acquired comfortably within one or two phases of the cardiac cycle would circumvent 
many of the diflRculties and residual artifacts associated with these compensation 
techniques. 

25 The speed with which magnetic resonance (MR) unages may be acquired has 

already increased dramatically over the past decade. The improvements in speed may be 
traced to a combination of advances in the technologies of magnet construction and 
actuation, and mnovations in imaging strategy. Strong, fast-switching magnetic field 
gradients and fast electronics have allowed the intervals between data collections to be 

30 reduced significantly. Meanwhile, fast gradient-echo and spin-echo sequences have 
reduced image acquisition time by allowing greater portions of A-space to be sampled 
after each spin excitation. Echo planar unaging (EPI), fast low-angle shot (FLASH), 
turbo spin echo (TSE), and spiral imaging techniques all allow very short intervals 
between acquisition of successive data points. The DUFIS, OUFIS, RUFIS, and 

35 BURST family of sequences fiirther reduce image acquisition time by eliminating time 
delays incurred during gmdient switching and echo formation. Details of the above- 
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mentioned eight techniques may be fomid in the following papers: P. Mansfield, Multi> 
planar image formation using NMR spin echoes . J. Phys. C. 10, L55-58 (1977); A. 
Haase, J. Frahm, D. Mattaei, W. Hanicke, K. D. Merboldt, FLASH imaging: rapid NMR 
imaging using low flip-angle pulses . J. Magn, Reson 67, 256-266 (1986); J. L. Listerud, 
5 S. Einstein, E. Cutwater, H. Y. Kressel, First principl es of fast spin echo. Magn. Reson. 
Q. 8, 199-244 (1992); C. Meyer, B. Hu, D. Nishimura, A. Macovski, Fast spiral 
coronary arterv imaging. Magn Reson, Med 28, 202-213 (1992); 1. J. Lowe, R. E. 
Wysong, DANTE ultrafast imaging sequence fPUFIS) . 1 Magn. Reson. Ser B 101, 
106-109 (1993); L. Zha, L J. Lowe, Optimized ultra-fast imaging sequence (OUFIS). 
10 Magn. Reson. Med 33, 377-395 (1995); D. P. Madio, L J. Lowe, Ultra-fast hnaging 
urin p low flip angles and FIDs. Magn. Reson. Med 34, 525-529 (1995); and J. Hennig, 
M. Hodapp, Burst imaging . MAGMA 1, 39-48 (1993). 

Increasing the speed of MR imaging further is a challengmg proposition, since 
the aforementioned fast imagmg techniques have already achieved an impressive 

15 eflRciency. All these techniques allow very short intervals between acquisition of 
successive data points, and hence do not waste much time in accumulating the data for 
the ifc-space matrix required to generate an image. In flow-encoded EPI images, for 
example, the entire complex A-space matrix is filled in a single spm excitation (which is 
followed by multiple spin conditioning cycles involving the application of multiple 

20 stepped field gradients), and the resulting im^e matrix is likewise "full " with useful 
information stored in both the real and the imaginary channels. One common feature of 
nearly all the fast imaging techniques currently in use, however, is that they all acquire 
data in a sequential fashion. Whether the requu^d data set, i.e., the A:-space data matrix, 
is filled m a rectangular raster pattern, a spiral pattern, a rapid series of line scans, or 

25 some other novel order, it is acquired one point and one line at a time. 

That is, the prior art in fast MR unaging has concentrated on increasmg the speed 
of sequential acquisition by reducing the intervals between scanned lines. Modifications 
to pulse sequences or to magnetic field gradients have yielded a gradual improvement in 
imaging speed by allowing faster sequential scanning of t-space, but these 
30 improvements face Ihnits due to the intervals necessary to create, switch or measure the 
nu^etic fields or signals mvolved in data acquisition. It would therefore appear 
difficuh to devise a sequential technique with significantly better efficiency than the 
current fast imaging techniques. 

Several fast imaging schemes have been proposed to date using simultaneous 
35 data acquisition in multiple RF coils, as described in: D. Kwiat, S. Einav, G. Navon, A 
decoupled coil detector arrav for fast image acquisition in magnetic res onance imaging. 
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Med Phys, 18:251-265 (1991); D. Kwiat, S. Einav, Preliminary experimental evaluation 
nf an inverse source imaging procedure using a de coupled coil detector array in 
ma gnetic resonance imaging . Med Eng Phys, 17, 257-263 (1995); J. W. Carlson, T. 
Minemura, Imaging time reduction through mult iple receiver coil data acquisition and 
5 ima ^e reconstruction . Magn Reson Med 29, 681-688 (1993) and U.S. Patent 4,857,846; 
and J. B. Ra, C. Y. Rim, Fast imaging using sube ncoding data sets from multiple 
detectors . Magn Reson Med 30, 142-145 (1993). These approaches have offered the 
promise of significant savings in image acquisition times. 

Carlson and Minemura describe a twofold acquisition time savings using 

10 two nested body coils. In their approach, partial data sets are collected 
simultaneously in the two coils, one of homogeneous sensitivity and the other 
with a linear gradient in sensitivity.Missing lines in t-space are generated usmg a 
series expansion in tenns of other phase-encoded lines. This approach using 
body coils appears to require that a significant portion of the data for the partial 

15 k-space matrix be acquired before any of the missing lines can be filled m by 
pos^rocessing, and thus does not allow for the missing lines to be built up as the 
data arrives, in real time. The approach uses coil sensitivity information in place 
of some portion of the gradient phase encoding steps, but has drawbacks. The 
coils used by Carlson and Minemura are body coils, which provide large volume 

20 coveiage but lower overall sensitivity than surface coils, and it would be difficult 
to augment their number to improve time savings. 

The approach of Ra and Rim involves a simultaneous acquisition technique in 
which images of reduced FOV are acquired in multiple coils of an array and the Nyquist 
aliasing in those images is undone by reference to component coil sensitivity 

25 mformation. The unaliasing procedure involves a pixel-by-pixel matrix inversion to 
regenerate the full FOV from multiple copies of the aliased image data. The 
"subencoding'* technique of Ra and Rim relies on estimates of component coil 
sensitivities by effectively probing the sensitivity at each pixel. This pixel-by-pixel 
approach can lead to local artifacts; for example, the matrix inversion can begin to fail in 

30 regions of low sensitivity. Further, by its very nature as a pixel by pixel dealiasing 
approach, the Ra & Rim method is computation-intensive and is limited to 
pos^rocessing, as all image data must be present before the reconstruction can be 
undertaken. 

In a related area, multiple coil signal collection has been used in MR phased array 
35 systems as reported m R.B. Roemer, W.A. Edelstein, C.E. Hayes, S.P. Souza, and O.M. 
Mueller, The NMR nhased array. Magn. Reson. Med 16, 192-225 (1990); C.E. Hayes 
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and P.B. P'^tw'^* - Nmse correlations in data simultaneousi v acquired from multiple 
>»..rface coil arrays , hdagn. Resort Med 16, 181-191 (1990); C.E. Hayes N. Hattes. and 
P.B. Pnftmftr, Volume imaging with MR phased arrays . Magn. Reson. Med. 18, 309-319 
(1991). The increased information content of the multiple receiyed signals in such 
5 systems has been used to increase the signal-to-noise ratio (SNR) of MR images. Since 
their initial description, phased arrays haye seen increasing use in clinical MR imaging. 
For example, the improyements in SNR proyided by phased arrays have allowed 
significant advances in imaging of the puhnonary vasculature as reported by T. K. F. 
Foo, J. R. MacFall, C. E. Hayes, H. D. Sostman, and B. E. Slayman, Pulmonary 
10 vasculature: si n gle hreath-hold MR imaein p with phased array coils. Radiology 183. 
473-477 (1992). Still, with a few notable exceptions, the bulk of phased array 
appUcations have addressed increased sensitfvity, with UtUe effort towards improving 
image acquisition speed or resolutioa 
Summary of the Invention 
15 An MRI system according to the present invention uses a multiple-coU data 

collection system to acquire some portion of the k-space matrix in parallel rather than 
sequentially in time. In a preferred embodiment, signals are obtained from multiple RF 
coils each occupying a different position with respect to the imaged volume, and each 
therefore having different but at least partiaUy overlapping spatial sensitivities. The 
20 multiple coils are preferably positioned and/or their outputs sampled in a manner to 
minimize inductive coupling, but they need not individually span the full region of 
interest nor be fully independent. The signals collected in this plurality of coils are then 
combined with suitably chosen weights to produce two or more composite signals, each 
of which approximates a sinusoidal sensitivity modulation, or "spatial harmonic." As 
25 used herein, tiie term "spatial harmonic" refers to a sinusoidal and/or cosmusoidal 
variation in spatial sensitivity witii a wavelengtii that is an integer fraction of tiie extoit 
of the field of view. Each line of these spatial harmonic composite signals constitutes an 
additional line of the A:-space matiix which would require a distinct gradient step in a 
conventional MR acquisition. Thus, tiie use of such signal combinations eliminates the 
30 need for some of the conventional gradient steps. The technique just outiined is referred 
to herein as SiMultaneous Acquisition of Spatial Harmonics (SMASH), and it may be 
used to reduce image acquisition times by a multiplicative factor witiiout a significant 
sacrifice in spatial resolution or signal-to-noise ratio (SNR). The invention also includes 
non-Fourier embodiments, in which tiie coil signals are tranformed or combined witii 
35 weights to yield composite signals which each correspond to a non-Fourier spatial 
distribution, such as a wavelet. 
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The SMASH technique in its preferred embodiment is implemented in an MR 
imaging apparatus and uses linear combinations of the simultaneously acquired signals 
fiom multiple surface coils to generate multiple data sets with distinct offsets in *-space. 
The full it-space matrix is reconstituted by interleaving these shifted data sets, and an 
6 image is therefore generated with only a fraction of the usual number of gradient phase 
encoding steps. Consequently, the total image acquisition time may be reduced or the 
total quantity of data generated during a fixed acquisition time may be increased by this 
same fraction. 

Conceptually, the invention may be seen as partially replacing gradient phase 
10 encoding by a spatial encoding procedure tied to the detection apparatus. In SMASH, 
some of the spatial modulations that normally distinguish differently phase-encoded 
lines are generated instead by amplitude modulation, resulting from combinations of 
component signals from multiple coils arrayed above or around die imaging volume. 
Computationally, by shifting responsibility for spatial encoding fiom tiie spin 
15 preparation stage to the stage of signal detection and combination, applicant is able to 
combine tiie multiple simultaneously acquired signals, eitiier immediately on-the-fly or 
after-tiie-fact, to produce multiple different spatial harmonics. In tiiis manner, multiple 
lines in Jfc-space are acquired simultaneously, in a parallel ratiier than a purely sequential 
acquisition scheme. 

20 A representative apparatus implementing flie invention simultaneously acquires 

partial signals fiom multiple coils in a surface coil array and combines tiiem into two or 
more differentiy weighted combinations tiiat accurately represent several spatial 
harmonics. The weights may be generated by theoretical calculations based on coil 
geometry, or else they may be derived by a calibration protocol carried out on a phantom 

25 image or on a pre-scan at die time of in vivo imaging. The calibration protocol uses a 
numerical optimization algoritiun to determine coefficients for linear signal combination 
which will best approximate the desired spatial harmonics. Once die spatial harmonic 
signal combinations have been formed, the missing lines in k-space arc filled in, and die 
conventional MR image is generated by Fourier tiransformation of tire expanded data 

30 matrix. Preferably, die coils m a SMASH coil array are surface coils, i.e. coils 

positioned on die surface of die body and designed to capture die MR signal efiBcientiy 
over a restricted region of interest Such surfece coils typically have a spatial extent on 
the order of 5-50 centimeters. 
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Brief Description of The Drawings 

These and other features of the invention will be understood from the description 
herein read in light of the prior art, together with the drawings illustrating representative 
embodhnents of the invention, wherein 

5 FIGURE 1 illustrates a general MRI apparatus having features common to the 

invention and to the prior art; 

FIGURE 2 illustrates a conventional data acquisition scheme; 
FIGURE 2A illustrates data acquisition of the present invention 
FIGURE 2B illustrates a receiver coil array of the invention and corresponding 
1 0 image plane geometry; 

FIGURE 3 illustrates another receiver coil array of the invention togetherwith 
coil sensitivities and the construction of spatial sensitivies from component coils; 

FIGURE 4 illustrates images formed from the coil array of FIGURE 2B and 
complex sensitivities ofthe coils; . 
1 5 FIGURE 4A illustrates fitting spatial hannonics with weighted coil sensitivity 

functions; 

FIGURES 5 and 5A show processing and representative images with the signal 
from a set of non-optimized coils and the fitting procedure of FIGURE 4A in a prototype 
embodiment of the invention; 
20 FIGURES 6A and 6B show MR images of a phantom made by processes ofthe 

prior artand the invention, respectively; and 

FIGURES 7A and 7B show MR images of a person's head made by processes of 
the prior artand the invention, respectively; 

Detailed Description 

25 FIGURE 1 illustrates schematically an MRI system 10 which includes 

the usual static magnet assembly, gradient coils, and transmit RF coils 
collectively denoted 12 under control of a processor 14, which typically 
communicates with an operator via a conventional keyboard/control workstation 
16. These devices generally employ a system of multiple processors for 

30 carrying out specialized timing and other functions in the MRI system 10 as will 
be appreciated. Accordingly, as depicted in FIGURE 1, an MRI unage processor 
1 8 receives digitized data representing RF NMR responses from an object region 
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under examination (e.g., a human body 1) and, typically vtamultiple Fourier 
transformation processes well-known in the art, calculates a digitized visual 
image (e.g., a two-dimensional array of picture elements or pixels, each of which 
may have different gradations of gray values or color values, or the like) which is 

5 then be conventionally displayed, or printed out, on a display 1 8a. 

In such overall operation, the apparatus of the present invention is 
largely conventional. However, in accordance with the present invention the 
basic RF data acquisition is modified, and subsequent signal processing altered, 
by providing a plurality of surface coils 20a, 20b...20i for simultaneous signal 

1 0 reception, along with corresponding signal processing and digitizing channels. 
The processor recombines the collected values into two or more spatial 
harmonics from which multiple lines of the signal matrix are developed. This 
recombination may be performed in real time as the data anrives, or after the feet 
via postprocessing as is convenient with the apparatus and the calibration 

1 5 information at hand. 

Operation of the invention will be better understood from a desoiption of 
gradient phase encoding in conventional processing and a comparison to the 
spatial encoding produced by amplitude modulation of the multiple coU signals 
in the {nesent invention, and a consideration of the information available in 

20 q>atial haimonics. 

In the genoal case, the ma gnetic re sonance signal for a plane with spin 
density lp(^>y)l and coil sensitivity may be written as follows: 



|s(k,.k,) = ffdxdyC(x,y)p(x.y)exp{-ik,x-ikyy}= C(H,.ky)® p(k„ky) 



25 where 



K -Y<^A| and \K ""y^y^l as usual, with g the gyromagnetic ratio, l^il 



and Eke magnitude of the x and gradients, and 0 and S the times spent the 
in thexaady gradients, respectively. Here, the spin excitation function as well 
as the effects of relaxation have been incorporated into a pulse-sequence-specific 
sensitivity fimction C. The symbol I®1 indicates a convolution. For regions of 
30 the sample in which the coil sensitivi ty is roug hly homogeneous we may write 



C(x,y)= 1| and |S(k^,ky)| [P^^'^^l the spatial Fourier transform of 



the spin density function. Double Fourier tr ansform ation with respect to El and 



0 reconstructs the usual spin-density image [ P^^y^ l 
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In a prior art phase-warp imaging method, a k-spaceinatrix of frequency- 
encoded and phase-encoded information, as represented goierally in FIGURE 2 
is produced for the image plane of interest, with successive lines represented by 
broad horizontal line corresponding to a different phase-encoding gradient value. 
5 Each phase encoding step corresponds to a distinct spatial modulation and these 
spatial modulations are represented schematically by the curves at the left of the 
FIGURE. 

Generally prior art magnetic resonance receiver coils, especially surface 
coils, do not have uniform sensitivity. Signals from different regions of the 
10 im^ed volume produce different currents in an RF coil, with the spatial 
variation in sensitivity being simply related to the inhomogeneity of RF field 
produced by the coil over the sample volume. For a standard circular suifece 
coil, there is a sensitivity "sweet spot" centered at roughly one diameter below 
the coil, with a monotonic falloff of sensitivity together with increasing phase 
15 differences in all duections. Traditional imaging protocols often position the 
receiving coil with the tai^et tissue at its region of maximum sensitivity. 

Applicant realized, however, that vAxen collecting a signal with some 
combination of coils having an aggregate sinusoidal and cosinusoidal spatial 
sensitivity profile, the MR signal fiom these coils will have an information 
20 content somewhat different fix>m that of the usual coil signal, and that by 
separating out one or more collected signals corresponding to pure spatial 
hannonics, these may be used to fill a larger portion of the data space than is 
conventionally done. A complex combination of cosinusoidal and sinusoidal 
sensitivity profiles yields 

C(x,y) = exp^Kyy}= cosKyy+isinKyy 



25 



[2] 



where iba is the spatial tKqaeacy of the inhomogeneous coil sensitivity. For the 
MR signal, we then have 



S(K<.^) = Jfdxdyp(x.y)exp{-ik^x- i(k, - K^)y}=p(k,,ls -Ky) 



[3] 



In other words, Ae combined MR signal from the inhomogeneous coils is 



30 shifted in A-space by an amount USJl This A-space shift has precisely the same 
form as the phase-encod ing shift produced by evolution in a y gradient of 
magnitude " Thus, when this shifted signal is acquired 
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simultaneously with the usual signal from a homogeneous coil, one gains an 
extra line in Jfc-space for each application of the x and y gradients. Applicant 
exploits this realization to convert the signals collected in a set of surface coils 
into spatial harmonic signals, and fill multiple lines of k-space from each signal 
5 collection. If one such harmonic is generated in addition to the homogeneous 
signal (zero^ harmonic) then one may then apply half as many y gradient steps 
to build up the full jt-space matrix necessary for image acquisition. When higher 
harmonics of larger 15J are added, the number of necessary y gradient steps is 
fiirther reduced. For a multiple coil setup representing M spatial harmonics m 
1 0 addition to the homogeneous signal, the acquisition time is reduced by an integer 
factor of If both positive and negative exponentials are used, 

C^(x,y) = exp^iKyy}= cosKyy± isinK^y 



[4] 



an extra line appears at an intervaJ^ above and below each gradient-encoded 
line. The *-space geometry of this acquisition strategy is illustrated in FIGURE 
15 2 A. 

Thus, if one were able to collect a signal with a coil shaped so that the 
response it picked up had the form of such a spatial harmonic, it would directly 
yield a k-space entry corresponding to one of the usual spin-conditioned 
measurements; an array of such coils might simultaneously acqmre several lines 
20 ofdata. 

Applicant further realized that the sensitivity of each individual coil in an 
array need not be strictly sinusoidal, so long as the net sensitivity of the array 
permits construction of harmonics havmg the desired sinusoidal shape. This 
realization significantiy relaxes the constraints on coil design and disposition, 
25 and in conjunction with a calibration or normalization process, allows flie signals 
from multiple coils with a wide range of shapes and geometries to be combined 
in several ways to yield net signals with multiple sinusoidal sensitivity profiles. 

FIGURE 2B illustrates geometry of the coil array and the image plane. A 
set of tiiree adjacent coils are used with separate coil ou^uts in a surface array 
30 extending generally over tiie region and planes of interest. The patient is 
positioned and tiie spin preparation fields are applied to condition a plane P 
which, by way of example may intersect tiie patient's heart to image the blood 
vessels tiiereof, or which may intercept a region of the abdomen to image its 
contents. The coil array 20 is located above or below the region to be imaged, so 
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that each coil 20a, 20b, 20c has at least some sensitivity to RF signals emanating 
from region P. In the illustrated embodiment the coils are each somewhat 
overlapped in the y-direction to minimize inductive coupling. As shown, they are 
overlapped at neighboring edges. Each coil, viewed individually has a sensitivity 
5 function which is highest directly above or below the center of the coil, and 
which falls off with distance from the coil center. In accordance with the present 
invention, the signals from these coils are combined to produce several separate 
"virtual" or synthetic signals which each correspond to a pure spatial harmonic. 
FIGURE 3 demonstrates this situation schematically for a set of eight 
1 0 rectangular coils 20a, 20b,...20h laid end-to-end, with a slight overlap. As shown 
in line (A) of the Figure, each coil 20a, 20b... has a sensitivity curve a, b,... which 
rises to a broad peak directly under the coil and drops off substantially beyond 
tiie coil perimeter. The sum of the coil sensitivities forms a relatively constant 
sensitivity, over tiie fiill lengtii of tiie array, corresponding to tiic zero* spatial 
15 harmonic. The remaining lines (B)-(E) of FIGURE 3 iUustrate lecombining of 
different ones of tiiese individual offset but otherwise identical coil sensitivity 
functions into a new syntiietic sinusoidal spatial sensitivity. Different 
weightings of the individual component coil sensitivities lead to net sensitivity 
profiles approximating several of the spatial harmonics of Ky. In tiie FIGURE, 
20 coil sensitivities (modeled schematically for tiiis Figure as Gaussian in shape) are 
combined t o produce har monics at various fractions of tiie fimdamental spatial 
wavelengtii K = ^"^^l witii E bemg on tiie order of tiie total coil array extent 



in>'. Weighted individual coil sensitivity profiles are depicted as tiiin solid lines 
beneatii each component coil. Dashed lines represent tiie sinusoidal or 
25 cosinusoidal weighting functions. Combined sensitivity profiles are indicted by 
tiiick solid lines. These combined profiles closely approximate ideal spatial 
harmonics across the array. 

For a weighted sum of component coil signals, the net sraisitivity profile 



C'°'(x,y)| is a linear combmation of tiie intrinsic sensitivity profiles of the 



30 component coils. In tiiis case, for an A^-component coil array, we have 

[5] 



C'^ix,y) = 'ZnjCjix,y) 



where nj is tiie weighting coefficient of tiie/* coil sensitivity fimction Cj(x,y). 
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In accordance with the aspect of Ae invention illustrated in FIGURE 3, 
by simply weighting the signals in the coil anray with appropriate weights and 
combining them, one obtains a signal whose amplitude is modulated by the 
spatial harmonic for that combination of the coil sensitivities. That is, the M (in 

5 this case, eight) component coil signals are acquired simultaneously and then 
multiplied by various weights and recombined in a total of up to N independent 
combinations each representing a spatial harmonic corresponding to a distinct 
gradient phase encoding step or distinct offset in k-space. This additional 
processing does not add in any way to image acquisition time, and it may be 

1 0 performed quickly on stored signal data as part of a post-processing algorithm or 
else may be implemented in a simultaneous flow-through fashion as described 
below. 

In actual practice, the sensitivity profiles of RF surface coils are not 
simple Gaussian profiles but more complicated fimctions which are, in general, 
15 complex in the mathematical sense of having both real and imaginary 
components. The coil sensitivity functions must describe both the magnitudes 
and the phases of the signals produced by precessing spins at various distances 
fiom the coil center, and these magnitudes and phases vary according to the 
reciprocity relation 



20 



W = -;|{B.^(r)«n(r)l 



[6] 



where E{r) is the voltage induced in a coil by a given voxel at position r, ni(r) is 



the nuclear magnetic moment of the voxel, and is the vector 

component of the field generated at r by a unit cunrent in the coU as described for 
example in D. I. Hoult and R. E. Richards, The signal-to-noise ratio of the 
25 nuclear magnetic resonance experiment J. Magn. Reson. 24, 71-85 (1976). 

This relationship has two main consequences for applicant's SMASH 
reconstruction: first, applicant assigns complex weights nj to capture the fiill 
information content of the coil s oisitivity functions, and second, the target 



complex exponential sensitivities 



C^{.x,y) = exp{±t^^} = GO&K^±ismK^ 



30 are preferably developed directly, rather than by separate generation and 
subsequent combination of sines and cosines. The real and imaginary parts of 
the complex coil sensitivity combina tions may be fitted simultaneously to the 



real and imaginary parts of F^^^'^^l- This may be done, for example, with a 
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numerical optimization algorithm using the sum of absolute magnitude 
deviations from the target profile or some other similar quantity as a measure of 
goodness of fit. In the prototype described below, tiiis fitting was performed m a 
straightforward manner and integrated vnth data reconstruction algontiuns 
5 written by applicant in tiie MATLAB programming language (Mathworks, 
Natick Massachusetts). Applicant has now confirmed tiiat such processing 
produces good image quality at greatly enhanced speed. 

In setting up an actual MR imaging apparatus, the coil geometry for 
SMASH signal acquisition may be that of an MR phased array, as described in P. 
10 B Roemer. W. A. Edelstein, C. E. Hayes. S. P. Souza, and O. M. Mueller, The 
MMi?nha.ed array . Magn. Reson. Med. 16, 192-225 (1990); C. E. Hayes and P. 
B pnrn-- ^nrrplritinns in dat a «mnlt;.nenuslv amnired fi^m multiple 

.nrf«..e. coil arrays. Magn. Reson. Med. 16, 181-191 (1990); and C. E. Hayes. N. 
Il.iUL-,-uiiir n Fnrmrr Vnlnir- —r-p MR pha.sed arrays. A/ggn. 
15 Reson. Med 18, 309-319 (1991). Indeed, many of the hardware components 
usefiil for SMASH imaging are aheady present in traditional phased arrays, 
which contain multiple inductively decoupled coils with some spatial separation, 
and which include separate receivers for independent collection of data from the 
coUs. The technical problem of minimizing inductive coupling of such coils has 
20 been previously addressed for multiple-coil constructions in tiie context of 

phased-array MR imaging devices, and two basic strategies were developed. One 
strategy is to design a coil array witii an appropriately chosen overiap of 
neighboring component coils to minimize inductive coupling. The second 
strategy involves the use of low input impedance preamplifiers on each 
25 component coU chamiel. Both of these features are advantageously applied in 
various embodiments of the prescaat invention. 
A Prototype and Results 

Applicant has implemented a basic embodiment of the invention as 
outiined above in tiie following manner. First, a linear RF coil array as shown in 

30 FIGURE 2B was selected having a geometry suitable for spatial harmomc 
generation. This array was a cardiac imaging array having tiiree eight-inch (200 
nmi) rectangular coUs adjacent to each otiier. The appropriateness of tiie coU 
geometry was first tested in numerical simulations, using analytic integration of 
tiie Biot-Savart law to calculate tiie transverse field Bx^r) of Eq. [6], and 

35 hence to model tiie sensitivity profile of each coil. Taking into account tiie 
number of spatial harmonics tiiat could comfortably be generated usmg tiie coil 
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array, and therefore the fraction of the total i-space matrix to-be collected, partial 
image acquisitions were planned, and a regimen of gradient steps and RF pulses 
were established. Image data were then acquired from each coil simultaneously 
in the separate channels of the array, using a conventional fast imaging sequence. 
5 Since only a fraction of the usual signal data was to be used to generate the fmal 
images, only a fraction of the usual imaging time was spent on data collection. 

Next, the stored individual coil image data were exported to a workstation 
for postprocessing. Individual complex sensitivity profiles were determined for 
each of the component coils in the array, using phantom image data instead of 
10 numerical simulations. This process is illustrated in FIGURE 4. For the 
phantom image technique a round water-filled dummy was suspended at the 
image plane used in the image acquisition above to provide a phantom image 
slice known to be of uniform spin density. Line (A) of FIGURE 4 shows the 
images a,b,c generated from the separate signal received in each of the respective 
15 coils 20a, 20b. 20c of FIGURE 2B, processed separately. Line (B) of FIGURE 4 
shows the corresponding real and imaginary components of the coil sensitivity 
function determined from the detected signal, taken along the diameter A of the 
FIGURE. 

To determine suitable weights for forming the desired spatial harmonic 

20 sensitivities from these three coils, an iterative fitting procedure was employed to 
represent each spatial harmonic as a linear combination of the three sensitivities. 
This was done as shown in FIGURE 4A. First, the real and unaginary 
components of the sensitivity fiinction were normalized, yielding the component 
values a,b,c plotted at panels A and B, respectively, of FIGURE 4A. Then, for 

25 the 200mm field of view, optimal weights for linear combination of component 
coil signals were determined by iterative fitting of the sensitivity data to the 
target spatial harmonic sensitivity profiles, with the complex weights serving as 
fitting parameters and the sum of absolute magnitude deviations from the target 
profile serving as a measure of goodness of fit Panel C of HGURE 4A shows 

30 the resulting best fit to the zero* harmonic, corresponding to uniform sensitivity. 
Panel D iUustrates the best fit to the first spatial harmonic. Each was realized as 
a simple linear combination, with complex coefiBcients, of the three coU 
sensitivities shown in panels A and B. 

In accordance with this aspect of the invention, the stored coil signals 

35 received during each spin conditioning cycle of image acquisition were 
recombined with the weights as determined above into two distinct spatial 



wo 9801600 



PCTAJS97/20068 



- 15- 



harmonics, conesponding to two distinct fc-space-shifted signals. The data sets 
were then interleaved, in the ordering sequence of FIGURE 2A , to yield the full 
ifc-space matrix. This matrix was double Fourier transformed in the conventional 
manner to give the reconstructed SMASH image. In the prototype, this data was 
5 reconstructed after the fact, i.e., post processed to both determine a suitable set of 
weights for fonning the composite spatial harmonic sensitivities, and then to 
lecombine the aheady measured signals with those weights into corresponding 
spatial harmonic measurement signals. In various embodiments of the device, 
however, the optimal weightings may be determined in advance withprior 
10 knowledge of the coil sensitivities, and the conversion of coil signals may be 
performed on-the-fly, that is, in real time. For example, a simple controller 
attached to the coils may determine the appropriate weights, from a lookup table 
or by some other linear or non-linear interpolative mechanism, and then combine 
the signals to produce a separate composite signal corresponding to each spatial 
15 harmonic. A neural netwoik may be used to "leam" the appropriate weighting. 

For simple coil configurations, the sensitivity may be modeled with 
accuracy, rather than be empiricaUy derived from normalizing the coil responses 
to a uniform sample. More generally, the process illustrated in FIGURES 4 and 
4A may be used to determine suitable coefficients for a given unage plane, and 
20 this process may be repeated to compile a table for each of many planes so that 
the apparatus need only look up the necessary weights in order to produce the 
desired composite signals for the plane of interest. When the sensitivity 
information is known in advance, then, reconstruction does not require a 
preliminary sensitivity measurement or an iterative fitting procedure, and simply 
25 involves a set of weighted sums. In either case, once the new composite signals 
are formed, the additional data must simply be entered at the correct matrbc 
position, and the fiill matrix subjected to a fast Fourier transform. This 
processing is as straightforward as the postprocessing which is commonly 
performed when conventional images are generated from phased arrays. 
30 Furthermore, once the fiill k-space matrix is constructed using the SMASH 
technique subsequent processmg of the k-space matrix to yield the image is 
identical to that presently employed in substantially aU Fourier MR imaging 
systems... 

FIGURES 5 and 5A illustrate the k-space matrix construction for this 
35 prototype in more detail. Half the usual number of gradient steps with twice the 
usual spacing in k-space were applied, and the RF response was recorded in each 
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of three coils. Given the gradient phase encoding regimen just described, each of 
these RF responses corresponded to an image of half the desired field of view. 
The coil signals were then combined into two sets with different weightings to 
produce synthetic composite signals corresponding to the zero* and first 
5 harmonics, and the two synthesized responses were interleaved to fill alternating 
lines of the k-space matrix. The matrix was then Fourier transformed in a 
conventional way to form the image. The resulting image had full resolution 
over the fiill field of view. As shown in FIGURE 5, line (A) depicts the 
acquisition of signal in each component coil. A representative signal point P is 
10 picked up in each coil, and thus contributes to the line in k-space formed with 
each coil signal. Line (B) shows thie processing performed in the prototype 
described above, according to which the coil signals arc combined into a first set 
So and a second set S, which each form different lines (even and odd as shown, 
corresponding to one gradient step separation) of the k-space matrix. These are 
1 5 combined at line (C) into the fiill matrix. 

FIGURE 5A illustrates corresponding images and signals. Line (A) of 
that FIGURE shows the half field-of-view images a, b, c reconstructed fiom 
separate coil signals, of which the signals themselves are shown in line (B). Line 
(C) shows the two spatial harmonic signals and line (D) the fiill interieaved 
20 signal, or k-space matrix. Line (E) illustrates the image reconstructed from the 
signals of line (C). The middle three stages (B), (C) and (D) correspond directly 
to the steps shown schematically in FIGURE 5. 

FIGURES 6A, 6B and 7A, 7B, respectively, show phantom and in vivo 
images formed by the three-coil prototype embodiment using the SMASH 
25 imaging technique. FIGURE 6A shows the proton-density image of a water- 
filled phantom, obtained using a TSE sequence (described fiirther below) with a 
total conventional acquisition time of 10 seconds. HGURE 6B shows the 
equivalent image obtained in 5 seconds with the SMASH processing of the 
present invention. Apart &om certain residual artifacts, the images have an 
30 identical appearance. 

Similar results were obtained in vivo. HGURES 7A and 7B show a 
reference and SMASH-reconstructed coronal image through the brain of a 
healthy adult volunteer, the images being acquired in 71 and 35 seconds, 
respectively. All these images were acquired using commercial hardware and a 
35 convenient pulse sequence with high spatial resolution and good SNR; a fifty 
percent acquisition time reduction was achieved. The acquisition time savings 
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described here have also been shown to apply to other conunercial machmes and 
pulse sequences, including the festest MR imaging machines and sequences. 
Some residual foldover artifects are present in the SMASH reconstructed image, 
due to imperfections in the composite spatial harmonic sensitivity profiles. 
5 However, the prototype used an existing coil array designed for other purposes, 
and improvements both in coil design and in the accuracy of coil sensitivity 
mapping and spatial harmonic generation, are expected to minimize these 
artifacts. 

The raw data for the foregoing images were generated on a commercial 
10 1.5 Tesla whole body clinical MR imager, which was a PhiUps NT, made by 
Philips Medical Systems of Best, Netherlands. For the phantom images, a 
standard circular phantom of 200 mm diameter and varied internal structure was 
used. A prototype coil array on hand in applicant's laboratory was used as the 
SMASH coil array. The equivalent c»cuit for this coil consisted of three 
1 5 rectangular component coils arranged in a linear array with a slight overlap in the 
right-left direction, similar to the array shown schematically m FIGURE 2B. As 
these coils were designed for another purpose, details of the coil circuitry not 
relevant here resulted in more complicated sensitivity profiles than would be 
expected from the simple geometry described above. It was these sensitivity 
20 variations which were primarily responsible for the residual foldover artifacts 
noted in FIGURES 6B and 7B. A simpler coil designed with spatial harmonic 
generation in mind will circumvent these artifacts. 

Further details of the prototype imaging process and equipment are given 

below. 

25 Image planning and data collection 

The phantom images of HGURE 6A, B were generated as follows. First 
with the phantom centered over the coil array, data for the reference image were 
acquired in a 6 mm thick coronal slice paraUel to and 80 mm above the plane of 
the array, using a turbo spin echo pulse sequence with five echoes per excitation. 

30 The field-of-view (FOV) was 200 mm, centered on the phantom, and matrix size 
was 256 x 256. Phase encodmg was performed in the right-left direction (i.e. in 
the direction of the coil array), with a single signal average. Data from each of 
the three component coil channels was acquired simultaneously and, in the 
prototype, were stored separately for later processing. Acquisition time was 

35 measured at 10 seconds. Next, a second coronal slice using the same technique 
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and imaging parameters was taken to serve as a measure- of component coil 
sensitivity. It was acquired 12 mm above the first one, in a region of uniform 
spin density in the phantom. Then, a third image at the same level as the first 
was obtained in half the time using twice the phase-encode step, and hence half 
5 the field of view, in the right-left direction. Matrix size was now 256 x 128. 
Acquisition time was 5 seconds, exactly half the time taken for the first image. 

For the in vivo images of FIGURE 7A, B, the volunteer was positioned 
with his head above the coil array, and images were taken in the same plane and 
ydth the same parameters as for the phantom images, except that eight signal 
10 averages with a slice thickness of 10 mm were used to improve SNR. The coil 
sensitivity image firom the phantom was also used as a sensitivity reference for 
the in vivo images. 

Determination of coil sensitivity and optimal component coil weights 

As desCTibed above m connection with FIGURE 4, coil weightings were 

15 determined by fitting actual coil sensitivity data to the desired spatial harmonic 
sensitivities. To do this, the sensitivity fimction along the right-left diameter of 
the phantom was extracted fi^om the component coil sensitivity reference images. 
The phantom spans the entire extent of each image, and since the sensitivity 
reference image plane intersects tiie phantom in a region of uniform spin density, 

20 these intensity profiles correspond precisely to the complex sensitivity fimctions 
of the coils along the diameter. The sensitivity of each component surface coil 
fells off monotonically with distance firom the coil, and the phase shift grows 
vwth distance horn the coil center. Using a gradient-descent fitting routine, these 
sensitivity fimctions were iteratively fitted to two target spatial harmonics: to the 

25 zeio-ftequency harmonic corresponding to uniform sensitivity, and to the first 
harmonic having a modulation wavelength equal to the field-of-view of 200 mm. 
The fit, shown in panels C and D of FIGURE 4A, demonstrates tiiat very close 
fits to tiie target harmonics may be achieved by component coil weighting, even 
witii only three coils. A slight residual ripple was present in die zero* harmonic 

30 fit and tiiis is primarily responsible for imperfections in the SMASH image 
reconstructions. 

Image reconstruction 

The reference images in FIGURES 6A AND 7A were formed by 
combining the component coil reference images pixel-by-pixel as tiie square root 
35 of the sum of square magMtudes. 
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FIGURE 5 A shows intermediate stages in reconstruction of the 

phantom image. The procedure for the in vivo image was identical. Using 
weights from the iterative fit, the three component coil signals (B) representing 
half-time, half-FOV aliased unages (A) were combined into two composite 
5 signal sets, one for the zero* spatial harmonic and one for the first spatial 
harmonic (C). Finally, the two composite signal data sets were interleaved to 
form a data matrix of size 256 x 256 (D), This matrix was Fourier transformed to 
yield the reconstructed image (E), which is also shown in FIGURE 6B. 

In general, the foregoing SMASH imaging technique may be applied with 

10 a great many of the known pulse sequences or spin conditioning techniques, and , 
will in general share the advantages of the underlying sequential imagmg 
methods which are used to collect partial ^-space information. As long as 
suitable spatial harmonics may be generated with a coil array, the additional 
acquisition time savmgs afforded by SMASH reconstruction involves no 

15 significant sacrifice in resolution or SlSfR. This contrasts markedly with the 
trade-off in SNR or resolution that characterizes many existing approaches to fast 
imaging, such as low flip-angle sequences. 

Signal-to-noise ratio 

Visual inspection of the phantom images in FIGURES 6A, 6B reveals a 

20 slight degradation in SNR in the SMASH reconstructed image as compared with 
the reference image. Part of this apparent loss may be traced to the residual 
foldover artifacts, in which some of the intensity of the primary image is "stolen'* 
by aliased ghosts. There is another noteworthy difference between the SNR 
profiles of the two images, however. The reference image in FIGURE 6A was 

25 generated using a sum-of-squares combination of component coil images, as is 
described in PB. Roemer, W.A. Edelstein, C.E. Hayes, S.P. Souza, and O.M. 
Mueller, The NMR phased array. Magn. Reson. Med. 16, 192-225 (1990) and as 
has become a standard practice in phased array imaging to improve SNR. This 
combination algorithm yields an essentially constant noise profile across the 

30 image, with varying signal (and hence varying SNR) that is enhanced in regions 
of significant overlap between component coil sensitivities. SMASH 
reconstruction, on the other hand, results in a constant-Jigwa/ unage with a 
varying noise profile. The linear combinations of component coil signals in 
SMASH are explicitly designed to produce a homogeneous composite signal. 

35 This is most obvious in the case of the zero^ harmonic combination, which 
produces a flat net sensitivity profile C^^^ across the image plane. While the 
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higher harmonics do involve significant spatial sensitivity variations, their 
profiles are all complex exponentials of unit modulus, and they do not lead to 
any intensity variations in the absolute magnitude image. (In other words, any 
apparent loss of sensitivity in the real channel is precisely compensated by a 
5 sensitivity gain in the imaginary channel.) SMASH reconstructions, therefore, 
do away with intensity peaks in regions of component coil overlap in favor of a 
spatially homogeneous image profile. Thus, sensitivity-dependent linear 
combinations may be used independentiy of SMASH as a metiiod of 
homogeneity correction. A homogeneity-corrected version of the image in 

10 FIGURE 6A may be generated by Unear combination of component coil 
reference images, using the weights calculated for the zero* harmonic profile in 
FIGURE 4A. Roemer et al. have described an alternative combination 
algoritiim for producing constant-signal images from phased arrays, and SMASH 
images produced with accurate spatial harmonics will have sunilar intensity 

15 profiles. 

Such considerations of innage intensity profiles are addressed primarily at 
the numerator of the signal-to-noise ratio. The behavior of noise m SMASH 
images is complicated by the presence of noise correlations among tiie variously 
combmed data sets that enter into a SMASH reconstruction. Noise fi-om each of 

20 the individual coils is repUcated (albeit with different weights) in each of tiie 
different harmonic composite signals, and tiiis might raise a concern tiiat some 
portion of die noise voltages will add coherently. If tije addition were purely 
coherent, i.e. if noise were to enter into each of the composite signals with 
identical weights, tiien tiie rms noise value would exceed die value for 

25 conventional full-time acquisitions by die square root of die number of 
composites. In otiier words, die upper limit of SNR loss due to noise correlations 
is the square root of the time savings fector. In practice, however, die signal 
combinations are formed widi distinct sets of complex weights, which reduces 
noise correlations substantially below this theoretical lunit 

30 Dependence on component coil sensitivities 

In die above described embodiment, die SMASH imaging technique 
depends upon die measurement and manipulation of coil sensitivities for its 
spatial hannonic generation. Errors in spatial harmonics may lead to aliasing 
artifects and SNR d^radation (for example due to mismatch and mixing between 

35 die spatial fiequency components of the image). This raises several practical 
issues. First, we must be able to estimate coil sensitivities with some reasonable 
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degree of accuracy. Sensitivity maps from phantom images were used for the 
SMASH reconstructions presented here. As noted above, numerical simulations 
of coil sensitivity may also be used, if these simulations have been reliably 
calibrated. The invention has also implemented use of other techniques for direct 
5 sensitivity estimation, including in vivo estimation using fast low resolution 
scans at the time of in vivo imaging. Low-pass-filtering of images could also be 
used to provide sensitivity maps. These techniques may be . automated for 
calibrating, normalizing and deriving weights for creating the composite signals. 

Another detail to be addressed is that Bq and Bj magnetic field 
10 inhomogeneities may distort the true coil sensitivity profiles and interfere with 
spatial harmonic generation. Spin echo sequences, like the TSE sequence used in 
this work, have the advantage of refocusing static field inhomogeneities. 
Gradient echo sequences have been found to perform just as welL The weight- 
finding algorithm used by applicant automatically compensated for some degree 
15 of residual field inhomogeneity and generated accurate spatial harmonics even in 
the presence of systematic phase or mtensity errors. In fact the technique has 
been found to be directly applicable not only to spin echo sequences such as TSE 
but also to several variants of gradient echo sequences including FISP and 
standard GRE. 

20 Finally, the geometry of the coil array will place certain limitations on the 

field of view (FOV), the position, and the angulation of image planes suitable for 
SMASH reconstruction. Surface coil sensitivity profiles vary with distance firom 
the coil center, and while simulations indicate that a wide range of image plane 
geometries are compatible with SMASH reconstruction, the reconstruction may 

25 begin to fail at large distances and angles, where sensitivity functions become 
broad and asynunetric. These constramts will be relaxed m the presence of larger 
numbers of component coils whose sensitivities provide good coverage of the 
imaging volume. The design of RF coil arrays specifically tailored for accurate 
and flexible spatial harmonic generation will help to overcome the limitations 

30 outlined above. 

In particular, the invention contemplates the use of coil arrays with 
multiple component coils extending m more than one linear direction. A two- 
dimensional array, such as an ^ xM rectangular array, will allow generation of 
spatial harmonics along multiple directions and will relax constraints on image 
35 plane position and angulation. Wrzqjaround arrays are also contemplated to allow 
spatial harmonic generation in a plane transverse tothe body plane. A coil 
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anangement with a surface array on the top of thew body "and another surface 
array on the bottom, possibly with some linear offeet with respect to one another, 
wUl allow fine tuning of spatial harmonics m a plane between the two and will 
have the added advantage of increasing overall SNR in such a plane. An 
5 extended grid coil grid coil may also be designed with a sensitivity profile more 
closely approximating a sinusoid. All such multiple-coil configurations will 
allow improvement in the generation of spatial harmonics, as well as allowing a 
greater number and variety of harmonics to be generated. As in the prototype 
embodiment described, the component coils in these arrays may be overlapped 
10 and may be ou^ut to low input impedance preamplifiers as necessary, to 
minimize inductive coupling. It is also possible that the simple expedient of 
transforming the coil sensitivities to achieve the spatial harmonic reconstruction 
procedure of applicant's SMASH invention may be applied in large anays of 
coilsof the sort proposed by MR Hutchinson and Raflf in 1988 and described 
1 5 further byKwiat, Einav and Navon, st^ra. 

It should be observed that, although recent studies have shown that the 
sensitivity advantages of phased arrays peak at coU numbers less than 10, the 
advantages of faster imaging with SMASH imaging may be further enhanced by 
arrays with larger numbers of coils and with new coil geometries. As one may 
20 mtuit from HGURES 3 and 4A, the higher the frequency of the spatial harmonic 
desired, the larger the number of component coils that will be required for its 
generation. This is in part because, at least for simple component coil shapes, the 
sharpest features of the net sensitivity fimction can only be on the order of the 
component coil dimension. As noted above.the extra sensitivity information 
25 from additional coils can also be used for the fine tuning of spatial harmonics and 
the improvement of image quality. 

The SMASH technique as applied to a conventional NMR apparatus 
described above partially replaces gradient phase encoding by a spatial encoding 
procedure tied to the detection coils, in which some of the spatial modulations 
30 that distinguish different phase-encoding lines are generated by combining 
signals from multiple coils arrayed above or around the imaging volume. This 
shift of responsibility from gradient geometry to coil geometry, and firom the spin 
preparation stage to the stage of signal detection and combination, allows for 
simultaneous acquisition of multiple lines of *-space. Enhancement of image 
35 speed, resolution or field of view by a fector of two, five, ten or more may be 
expected with appropriate coil anrays generating sufficient harmonics. 
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Apart from its promise for progressively faster acquisitions, SMASH has 
a number of practical advantages as a fest imaging scheme. As a partially 
parallel acquisition strategy, it can be combined with most existing sequential 
fest imaging techniques for multiplicative time savings. No special hardware is 
5 required, other than an q)propriate coil array. 

One embodiment of the invention includes a coU array together vwtii a 
digital signal processor configured to , combine the outputs of the component 
coils with appropriate weights as described above, and produce two or more 
output signals, each of which represents a composite spatial hamionic as 
1 0 described above. The array vwtii processor may then be directiy substituted for a 
conventional receiving coil, with the difference tiiat it produces two or more lines 
of data for each spin preparation cycle of the MRI apparatus. As such, the plug- 
in coil unit may operate with older MRI devices to directiy enhance the 
acquisition time, field of view, or resolution by a fector of two or more, without 
15 changing the expensive magnet and otiier spin preparation hardware of the 
device. CommerciaUy available or home-buih MR phased arrays may well 
suffice to yield significant time savings, and SMASH may be performed on 
machines not equipped with much more costiy magnetic enhancements such as 
EPI gradient systems. Combination of tiie component coil signals may be 
20 performed after the fact, allowing for a wide range of postprocessing steps, 
including fine tuning of spatial harmonics, adaptive artifact correction, or even 
non-Fourier (e.g. wavelet) encoding and reconstruction. Furthermore, it should 
be noted that the invention fiirther contemplates that both for constructing tiie 
basic component coil transforms or weights and for any tuning or adaptive 
25 adjustments, tiie invention may include a digital signal processor or neural 
network processor in addition to the normal complement of control and 
processing assemblies to evaluate the responses or coU signals and cany out the 
signal combination into spatial harmonics in real time. Even in tiie simplest 
embodiments, spatial harmonic reconstructions automatically yield a 
30 homogeneous intensity profile for tiie reconstructed image, and tiiis may be 
advantageous for some imaging j^jplications. 

The time gained in a SMASH acquisition may be used to collect extra 
data for better spatial resolution and/or SNR, or else faster acquisitions may be 
used to eliminate motion artifects fi»m mobile structures in tiie field of view 
35 The invention has been described above witii reference to a prototype 

embodiment which produces signals corresponding to spatial harmonics tiiat 
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yield additional k-space lines equivalent to the conventioifiar gradient steps in a 
Fourier transfoim imaging system and a simple coil airay and processing 
assembly for enhancing image speed, quality or field. However, the method is 
broadly applicable to coils and processors for producing various spatial 

5 representations, including wavelets or other non-Fourier spatial representations to 
replace, supplement or enhance MR imaging processes and systems. The 
invention being thus described, further variations and modifications will occur to 
those skilled in the art, and all such variations and modifications are considered 
to be within tiie scope of die invention and its equivalents, as set forth in the 

1 0 claims appended hereto. 



What is claimed is: 
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Claims 

1. A method of imaging a region of a body by MR by capturing MRl data during 
multiple cycles of NMR spin conditioning and measurement, said method comprising 

5 the steps of: 

i) during one spin conditioning cycle receiving an NMR RF measurement 
signal response simultaneously in each of plural surface coils having different spatial 
sensitivities, 

ii) combining said NMR RF measurement signal responses received in the 
1 0 plural coils into plural sets of composite NMR signals, each set of composite NMR 

signals containing different spatial information, 

iii) repeating said steps of receiving and combining during additional spin 
conditioning cycles to form a data set augmented by at least one set of said composite 

NMR signals, and 
15 iv) foiming an MR image firom said data set 

2. A method according to claim 1 , wherein said steps of receiving and combining 
are carried out to form images of a given field of view and spatial resolution while 
reducing at least one of number or duration of spin conditioning cycles and thereby 

20 reducing image acquisition time. 

3. A method according to claim 1, whercm said steps of receiving and combining 
are carried out to form images in which at least one of, field of view, spatial resohition, 
and acquisition speed is oihanced. 



25 



4. A method according to claim 1, wherein tfie step of combining combines the 
measurement signal responses in linear combinations. 



5. A metiiod according to claim 4, wherein the Unear combinations are performed 
30 using complex coefficients. 

6. A method according to claim 1 , claim 4, or claim 5, wherein tiie step of 
combining signal responses yields composite spatial sensitivities which approximate 
spatial harmonics of an imaged field of view. 

35 



« 
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7. A method according to claim 6, in which the spatial termonic composite signals 
are substituted for data lines corresponding to omitted increments in magnetic field 
gradients. 

5 8. A method according to claim 6. wherein the step ofreceiving in surface coils 
includes receiving in at least one of a two dimensional coU array, mNxMcoil array, a 
wrap around coil array and an extended grid coil array. 

9. A method according to claim 1 , further comprising the step of sensitivity 

1 0 mapping said plural coils and determining a transform for combining the received NMR 
RF measurement signal responses into said sets of composite NMR signals. 

10. A method according to claim 6. wherein the cycles of spin conditioning and 
measurement are selected firom among spin echo, gradient echo and BURST style 

15 ima^g techniques. 

11. A method according to claim 1 , wherein the step of combining into sets of 
composite NMR signal responses combines into sets of signals corresponding to non- 
Fourier spatial distributions, and the step of processing applies a corresponding non- 
20 Fourier transformation to form the MR image. 

12. A method according to claim 1 , wherein the step of combining into sets of 
composite NMR signal responses combines into wavelets forming sets of spatiaUy 
orthogonal data, and the step of processing applies a wavelet transform to form the 

25 image. 

13. A magnetic resonance imaging (MR!) apparatus of the type having a magnet for 
establishing a magnetic field in a region means for conditioning nuclear spins in the 
region and means for collecting an RF signal fi-om the conditioned spins, said magnet 

30 and said means being operated in repetitive cycles to collect spin data and construct an 
image of material in said region, wherein said means for coUecting an RF signal includes 
a plurality of surface coils forming an array, and each having a different 

spatial sensitivity, and 

means for combining signals firom plural individual coils of said array 
35 with different weights into a plurality of sets of signals corresponding to orthogonal 
^tial representations of spin data firom said region. 
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1 4. A magnetic resonance imaging (MRI) apparatus according to claim 13, wherein 
said means for combining combines said signals in jflow-through fashion as said spin 
data are collected.. 

5 15. A magnetic resonance imaging (MRI) apparatus according to claim 1 3, wherein 
said spatial representations are spatial harmonics. 

16. A magnetic resonance imaging (MRI) apparatus according to claim 1 3, wherein 
said spatial representations include a non-Fourier representation. 

10 

1 7. A magnetic resonance imaging (MRI) apparatus according to claim 1 6, wherein 
said non-Fourier representations are wavelets. 

1 8. A magnetic resonance imaging (MRI) apparatus according to claim 13, further 

1 5 comprising means for determining said weights, and wherein said means for combining 
combines in accordance with weights determined by said means for determining. 

19. A magnetic resonance imaging (MRI) apparatus according to claim 1 8, wherein 
said means for determining includes a stored table of weights for each of a plurality of 

20 imaging regions for forming said spatial representations. 

20. A magnetic resonance imaging (MRI) apparatus according to claim 13, wherein 
said plurality of surface coils form an array having an extent of approximately twenty to 
about fifty centimeters, and individual coils of the array have an extent each have an 

25 extent under about twenty centimeters. 

21 . A magnetic resonance imaging (MRI) apparatus according to claim 1 3, wherein 
said plurality of surface coils form a linear array. 

30 22. In a magnetic resonance imaging (MRI) apparatus of the type having a magnet 
for establishing a magnetic field in a region, means for conditioning nuclear spins in the 
region and means for collecting an RF signal firom the conditioned spins, said magnet 
and said means being operated in repetitive cycles to collect spin data and construct an 
image of material in said region, the improvement wherein said means for collecting an 

35 RF signal includes 

a plurality of surface coils forming an array, and each having a different 

spatial sensitivity, and 
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means for combining signals from plural individual coils of said array 
with diiferent weights into a plurality of sets of signals corresponding to orthogonal 
spatial representations of spin data from said region. 
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frequency encode (Kx) 

FIGURE 2 cJepIcllon of Ihe standard procedure for MR phase-warj 

imaging. Each horizontal line in the signal data set con-esponds 
to a different phase-encoding gradient value. The phase-encodir 
gradients produce spatial modulations shown as Insets on the left 
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Iterative fit of spatial harmonics with weighted coil sensitivity functions. 
A,B: Individual coil sensitivity functions nomialized to unit absolute 
magnitude. C: Best fit to homogeneous sens!ti>4^. (zeroth hannonic). 
D: Best fit to first harmonic. 
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